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Abstract 



A theoretical investigation of quantum-transport phenomena in meso- 
scopic systems is presented. In particular, a generalization to "open systems" 
of the well-known semiconductor Bloch equations is proposed. The presence 
of spatial boundary conditions manifest itself through self-energy corrections 
and additional source terms in the kinetic equations, whose form is suitable for 
a solution via a generalized Monte Carlo simulation. The proposed approach 
is applied to the study of quantum-transport phenomena in double-barrier 
structures as well as in superlattices, showing a strong interplay between 
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phase coherence and relaxation. 
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The Monte Carlo (MC) method, which has been applied for more than 25 years for 
calculation of semiclassical charge transport in semiconductors, is the most powerful nu- 
merical tool for microelectronics device simulation However, present-day technology 
pushes device dimensions toward limits where the traditional semiclassical transport theory 
can no longer be applied, and a more rigorous quantum transport theory is required [Q. 
To this end, various quantum-kinetic formulations of charge transport have been proposed, 
based on Green's function |§ or Wigner-function f| approaches. While such quantum- 
mechanical formalisms allow for a rigorous treatment of phase coherence, they typically 
describe energy-relaxation and dephasing processes via purely phenomenological models. A 
full quantum-mechanical simulation scheme for the analysis of transient-transport phenom- 
ena in the presence of carrier-phonon interaction has also been proposed |J. However, due 
to the huge amount of computation required, its applicability is still limited to short time- 
scales and extremely simplified situations. As a result, despite many efforts and despite 
the unquestionable intellectual progress associated with the study of these quantum-kinetic 
formulations, their application to realistic devices in the presence of a strong scattering dy- 
namics is still an open problem. Recent results by Datta, Lake, and co-workers seem to be 
rather promising |5[]. However, their steady-state Green's function formulation cannot be 
applied to the analysis of time-dependent non equilibrium phenomena, which play a crucial 
role in modern optoelectronic devices. 

In this letter we propose a generalized MC approach for the analysis of hot-carrier trans- 
port and relaxation phenomena in quantum devices. The method is based on a MC solution 
of the set of kinetic equations governing the time evolution of the single-particle density 
matrix. Our approach can be regarded as an extension to open systems of the generalized 
MC method recently proposed for the analysis of the coupled coherent and incoherent carrier 
dynamics in photoexcited semiconductors |7|] . Compared to more academic quantum- kinetic 
approaches []3|-[5] — whose application is often limited to highly simplified physical models 
and conditions — , the proposed simulation scheme allows to maintain all the well known 
advantages of the MC method in describing a large variety of scattering mechanisms on a 
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microscopic level Q. 

In order to properly describe carrier-transport phenomena in mesoscopic structures, an 
electron-phonon system can be considered, whose Hamiltonian can be schematically written 
as H = H + H'. Here, the single-particle term H D includes the phonon and free-carrier 
Hamiltonians as well as the potential profile (including possible external fields), while the 
many-body contribution H' accounts for all possible interaction mechanisms, e.g. carrier- 
carrier and carrier-phonon coupling. By denoting with a (r) = (r|a) the wavefunctions of 
the single-particle states a |§ and with e Q the corresponding energies, the equation of motion 
for the single-particle density matrix p |§ in this a-representation can be schematically 
written as: 



d d 



Ho dt 

The time evolution induced by the single-particle Hamiltonian H Q can be evaluated exactly. 
In contrast, the contribution due to the many-body Hamiltonian H' involves phonon-assisted 
as well as higher-order density-matrices; thus in order to "close" our set of equations ap- 
proximations are needed. In particular, as described in |TTJ, the "mean-field" approximation 



together with the Markov limit leads to a closed set of equations still local in time. Within 
such approximation scheme, the equations of motion in (|]) can be written as 

"JlPoiP = ^2 Lap^a'P'Pa'P' (2) 
at a ,p, 

with 

L a /3,a'l3> = TZ ( e a ~ Cs) ^a/3,o'/3' + ^a/3,a'/3' ■ (3) 

Here, the two terms correspond to the separation in (|1|). The explicit form of the scattering 
tensor T involves the microscopic in- and out-scattering rates for the various interaction 



mechanisms ||11|| . The above kinetic equations can be regarded as a multiband generalization 
of the well known semiconductor Bloch equations (SBE) ||12|| . 

The analysis presented so far is typical of a so-called "closed" system, i.e. a physical 
system defined over the whole coordinate space. However, this is not the case of interest 



for the study of quantum-transport phenomena in mesoscopic devices, where the properties 
of the carrier subsystem are strongly influenced by the spatial boundaries with the external 
environment. This requires a real-space description, which can be obtained in terms of 
the phase-space formulation of quantum mechanics originally proposed by Wigner [13|] and 



generalized to solids in the pioneering paper by Buot ]T3j]. In our case, this corresponds to 
introduce the following unitary transformation u connecting our aft representation to the 
desired phase-space r, k: 

u aP (v, k) = (27r)-§ J dT><P a (r + ±r') e-* W ^ (r - ~r') . (4) 

By applying the above Weyl transform to the single-particle density matrix p, we obtain the 
so-called Wigner function [Q: 

f w (r, k) = J2 Papu a p(r, k) . (5) 

a/3 

For a closed system, the Wigner function f w is defined for any value of the real-space 
coordinate r and its time evolution is fully determined by its initial condition. In contrast, 
for an open system f w is defined only within a given region Q of interest and its time 
evolution is determined by the initial condition inside such region as well as by its values 
f b w on the boundary of the domain at any time t' > t Q . More specifically, by applying the 
Green's function theory to the equation of motion for f w — which is obtained by applying 
to Eq. (|2|) the Weyl- Wigner transform ([5]) — we get: 

f w (r, k; t) = jf dr' J dk!G{v, k; r', k'; t - t Q )f w (r', k'; t ) 

+ Jdr b J dk' jf* dt'G(r, k; r 6 , k'; t - t')f^(r b , k' , t')v(k') , (6) 

where 

G(r,k;r',k';r)= £ u a0 (v, k) [e L A „ ,„,<^(r', k') (7) 

is the evolution operator, while v{k) is the component of the carrier group velocity normal to 
the boundary surface. We clearly see that the value of f w is obtained from the propagation 



of the initial condition f (t ) inside the domain Q plus the propagation of the boundary 
values f b w from the points of the surface at any time t' to the point r, k of interest. 

Given the above Wigner formulation for open systems, we now introduce a corresponding 
density-matrix description via the following "inverse" Weyl- Wigner transform |T5| : 



Pa/3 



Jjr J dku* a p(r,k)f w (r,k) . (8) 



By applying the above transformation to Eq. (^) and then performing its time derivative, 
we finally obtain: 

= L aP,a>l3>P a ipi + S af3 , (9) 

al a'P' 

where L = ULU~ l is the Liouville tensor @ "dressed" by the transformation 

U a p, a 'i3' = J n dr J dk<^(r,k)w Q ,/3'(r,k) , (10) 

while 

S af3 = u af3,a>f3> / dr b j dku* a , p ,{v b , k)v(k) (r b , k) (11) 

a'P' J J 

is a source term induced by our spatial boundary conditions. 

Equation (|J) is the desired generalization to the case of open systems of the standard 
SBE in Eq. (|j). In addition to the source term in Eq. (|TT|) , the presence of spatial boundary 
conditions induces modifications on the Liouville operator L of the system via the transfor- 
mation U in Eq. (|T0|). 

The generalized SBE @ can be still regarded as the result of a coherent single-particle 
dynamics plus incoherent many-body contributions [see Eq. ([I])]. Therefore, they can be 
solved by means of the same MC simulation scheme described in [0] . The method is based 
on a time-step separation between coherent and incoherent dynamics: The former accounts 
in a rigorous way for all quantum phenomena induced by the potential profile of the device 
as well as for the proper boundary conditions. The latter, described within the basis given 
by the eigenstates a of the potential profile [IH], accounts for all the relevant scattering 



mechanisms by means of a conventional "ensemble" MC simulation 0. 
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In order to illustrate the power and flexibility of the proposed theoretical approach, we 
have simulated quantum-transport phenomena in rather different physical systems, namely 
double-barrier structures and superlattices. Since we are interested in low temperature 
and low carrier density conditions, only optical-phonon scattering has been considered. We 
have first carried out the simulation of an electron wavepacket entering the double-barrier 
structure [17| of a GaAs/AlGaAs resonant tunneling diode (RTD). Figure |l] shows the time 



evolution of the wavepacket in the absence of scattering as a function of position (a) and 
energy (b). It is easy to recognize the well-established resonance scenario typical of any 
purely coherent dynamics: as the wavepacket enters our RTD structure, a part of it is 
transmitted and a part is reflected [see Fig. |l|(a)]. Since in this simulation scattering is 
not included, the wavepacket central energy is conserved, i.e. no energy relaxation occurs 
[see Fig. 0(b)]. On the contrary, in the presence of incoherent scattering processes the 
resonance dynamics of Fig. |I|(a) is strongly modified by the scattering itself, as shown in 
Fig. 0(a). In particular, the presence of phase-breaking scattering processes is found to 
reduce both the interference peaks and the transmitted wavepacket. This is confirmed by 
the corresponding energy distribution in Fig. 0(b), where we clearly recognize the granular 
nature of the dissipation process through the formation of optical-phonon replica. This 
is the fingerprint of any full microscopic treatment of energy relaxation, thus confirming 
the microscopic nature of our quantum-mechanical simulation, in contrast with all previous 
phenomenological approaches. 

As a second testbed for our generalized MC approach, we have also simulated electrically 
injected Bloch oscillations in superlattices (SLs). The system under investigation consists 
of a biased GaAs/AlGaAs SL [18[ surrounded by two semi-infinite GaAs regions. In our 
simulated experiment an electron wavepacket is injected from the left contact (GaAs region) 
into the SL region [see Fig. [| . Figures [3] and £| illustrate the time evolution of the wavepacket 
with and without scattering, respectively. When the wavepacket reaches the SL structure 
most of it gets reflected backwards, while some portion of it tunnels into the Wannier-Stark 
ladder associated with the SL and starts to oscillate at a frequency of about 3.5 THz. Every 



time the packet reaches the boundary of the oscillation region a part of it gets transmitted 
via Zener tunneling. We should notice, however, that such Zener processes do not destroy 
the Bloch-oscillation dynamics, but simply reduce the charge density within the SL region. 
Indeed, in the scattering-free case [see Fig. |3[] the Bloch oscillations are found to persist 
on a picosecond time-scale. In contrast, once scattering mechanisms are considered [see 
Fig. U, the phonon- induced dephasing drastically reduces their lifetime. As for the case of 
the RTD, the carrier dynamics is the result of a strong interplay between phase coherence 
and relaxation. 

In conclusion, we have proposed a generalization to open systems of the well-known SBE. 
This approach allows for a proper description of the strong coupling between coherent and 
incoherent dynamics. Indeed, our simulated experiments clearly show the failure of any 
purely coherent or incoherent approach in describing typical quantum-transport phenomena 
in semiconductor nanostructures. 

Contrary to all previous quantum-transport investigations, the proposed theoretical 
scheme allows to fully recover — and benefit from — the unquestionable advantages of the 
semiclassical Monte Carlo simulation, thus opening the way to the theoretical modelling of 
realistic quantum devices. 

We are grateful to Carlo Jacoboni, Tilmann Kuhn, Roger Lake, and Massimo Fischetti for 
stimulating and fruitful discussions. This work was supported in part by the EC Commission 
through the TMR Network "ULTRAFAST" . 
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FIGURES 

FIG. 1. Carrier density at different times as a function of position (a) and energy (b) cor- 
responding to an electron wavepacket injected into a RTD structure in the absence of scattering 
processes (the two barriers are schematically depicted as black vertical lines). 

FIG. 2. Same as in Fig. [I], but in the presence of scattering processes (see text). 

FIG. 3. Contour plot of the charge density corresponding to a wavepacket electrically injected 
into a finite SL region (marked with vertical white lines) in the absence of scattering processes. 

FIG. 4. Same as in Fig. H, but in the presence of scattering processes (see text). 
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